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The collective description of the octupole degrees of freedom has been a long standing problem in nuclear physics [0 . 
The phenomenological description of low-energy vibrational states |^,|| discussed them in terms of simple "surface 
modes" thus allowing for the description of the octupole degrees of freedom in terms of seven collective variables 
asfj, [^,||. In the eighties, theoretical calculations |5| predicted the existence of octupole stable deformations and 
1 aroused considerable interest, especially in the Ce-Ba and the Rn-Th regions. The level scheme of a 
. few moderately- or weakly- deformed nuclei, such as ^^Ge ^^^Sm |Q,^^®i?a ||^ or ^"^^'^^^iJa |9| presents features 
T-H that may be related to octupole instabilities and softness of the nucleus with respect to possible exotic octupole 
■ deformations. Lately |lO| there has been evidence for the existence of stable octupole deformations in the Rn-Th 
region making the problem of the collective description of this degree of freedom more actual. It has been shown that 
this type of collective excitations can significantly change the fusion cross section for heavy ions Initially, only 
octupole deformations of the ¥30 type were considered, but the possibility of other type of deformations (such as Y32 
Ip^or F31 has also been discussed. The study of more general structures was intended without much success @. 
There were also several proposals related to the appropriate parametrization of the collective variables. Rohozinski 
^ ' proposed a parametrization based on the symmetries coming from the Oh group of transformations of the frame of 

(~| , reference. Hamamoto [ p^ wrote octupole shapes in terms of the irreducible representations of the octahedron group 
^ Ai, Fi{k) and F2{k) and suggested that a parametrization that covers all possible octupole deformed shapes (without 
. ^ , double counting the same shape) is provided by removing F2{k) (eliminating in this way three parameters). 
^ ■ In the present paper we return to the idea used for the definition of the intrinsic system in Bohr's paper |^ (also 
5h [ used by other authors for the collective description of the pairing degrees of freedom): in the intrinsic system the 
. . 1 inertia tensor (related to the rotational part of the kinetic energy) must be diagonal. 

In Sect. 2 we review some features of the collective excitations, the Hamiltonian, we discuss how to separate the 
rotational and intrinsic degrees of freedom and we propose a simple parametrization. We use the Pauli prescription to 
quantize the quadratic Hamiltonian. In Sect. 3 we discuss the structure of the wave functions and of the elementary 
scalars or invariants (as defined in Rcf. [jl^), as well as the symmetries of the wave function. In Sect. 4 we discuss 
the rigid rotor that illustrates the differences between the rotational bands related with octupole deformations as 
compared with the quadrupole ones 118]. Section 5 summarizes the results obtained in the present paper. 
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II. FORMULATION OF THE PROBLEM 



A. The collective variables 



The theory of collective oscillations has been developed long time ago by several authors H,^,^ . The classical case 
corresponds to quadrupole oscillations that were described in detail by A. Bohr in 1952. A very clear and pedagogical 
discussion concerning the symmetries and wave functions for this problem can be found in a paper (l| by K. Kumar 
and M. Baranger, where a numerical solution for this problem was first proposed. A collective treatment was also 
performed for the pairing degrees of freedom: the pairing acting only on one type of particles yields a two dimensional 
collective description ]20| |, while the corresponding to the T = 1 case yields a 6 dimensional one [ plf . In the three 
cases described, the recipe used to found the collective description has the same basic ingredients: 

a) The kinetic energy was written explicitly in terms of the time derivatives of the collective variables. 

b) The collective variables (of dimension five, two or six, depending on the case) were written in the intrinsic system, 
isolating the intrinsic variables from those that describe the motion of the intrinsic system (the three Euler angles in 
the five and six dimensional cases, one angle in the two dimensional one). In all cases the intrinsic system was defined 
as the system where the inertia tensor (the part related with collective rotations) were diagonal. 

c) The kinetic energy was written in terms of the Euler angles and the intrinsic variables, and the coupling between 
both types of degrees of freedom was studied. Regarding the potential energy it has parts that depend only on 
the collective variables (and not on their time derivatives). This part can be written in terms of elementary scalars 
p7| that are constructed using the collective variables. For example, in the quadrupole case there are only two of 
such scalars (/3^ and cos (87), using Ref. [D notation), in the pairing case (between particles of one type) there is 
only one of such scalars (A^ using Ref. pO| notation) and in the T = 1 pairing case there are three of such scalars 
(A^, A2e±'**'*cos(2r), using Ref. Q notation). 

d) The general symmetry properties of the wave functions that are related to different choices of collective (Euler) 
angles and intrinsic variables corresponding to the same collective variables in the Lab. system are used to define 
the range of the different intrinsic variables. In some simple cases it was also possible to write down explicitly some 
elementary tensors 22 2^ in terms of the collective variables. 



In the present paper we will follow this program for the octupole surface oscillations, relating them to the shape of 
the nucleus. Many times the surface of the nucleus, in polar coordinates, has been expressed in terms of the nuclear 
radius as 

i?(0, 0) = Roil + «A,M>A,^(0, 0)) , (2.1) 

where Rq is the radius of the nucleus in its spherical equilibrium shape and cka,/^ are the collective coordinates that 
describe the deformation of the nuclear surface. We use for the Y\^^ the spherical harmonics satisfying the Condon and 
Shortley convention, i.e. ^a.-z^ = (—1)^ (Yx^p,)* . As the radius has to be real, it follows that 0;^,-^ = (—1)'^ (q^a.^i)*- 

It is assumed that these variables change slowly with time, and therefore it is usual to express their kinetic energy 
as a quadratic function of the velocities as 
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(2.2) 



where is the time derivative of a^. It is well known that the coefficients ax.fi written in Lab. system and the ax.fi 
written in the intrinsic system (that we will denote by ax.fi to avoid confusions) are related by (in what follows we 
use the notation of Ref. p9[|) 



ax 



where cui are the three Euler angles. As usual the I?^ ^(<^i) functions, which are related to the matrix element of the 
rotation operation, are defined by 



T^lA^^) =< A,Af|exp(ii/'/2)exp(i6i/y)exp(i(/)4)|A,i^ >, (2.4) 



where uji = {i>, 9 , (j)) are the Euler angles. 



From Eq.(|2.3D it follows that 



2 



a>,,^ = Y^[D^^^,{oJi)]*ax,., (2.5) 

with = 

In order to evaluate the kinetic energy we must write down exphcitly the time derivative of ct\,n- 

= Y.KA^^)T<- + YPi'^^^^)^*^^^- ■ (2-6) 

V If 

The time derivative of the functions can be expressed as 

^lA^i) =iJ2l-'^irn{>^i^Wu,m , (2-7) 

where -M„ are the (2A + 1) dimensional representation of the angular momentum operator and satisfy the commutation 
relations 

[Mi,M2] = -iM3, (2.8) 

and cyclic permutations (note the signs). It is customary to use a representation where (M3)m,n = mdm^n- The 
magnitudes are the angular velocity components in the k direction. We can then express the kinetic energy as 

^ 

^ ~ ^ ^ (Q^A,//) ~ ^vibration ~i~ ^rotation ~i~ ^coupling ? (2-9) 



where for simplicity a single A is considered, and 

B 



fi,i^.v .m,m ,K,K 

Tcoupiing=ij Yl [-«A,.'«A,.2?^:„^2?^,,'9«(M,),,„ + aA,.a*_^,P^;,I?^,„g4M,)„,,]. 

Using the unitarity condition for the 'D'^^^{uji) functions ■^Ai*i'('^i)^^ ('^i) = K,v'\^ the vibrational part of the 
kinetic energy simplifies to 

'^vibration — ~ ^\,v^\,v i (^-H) 



while the rotational part is 
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%-otaUon = y X) Y1 °'>'^''°'ly (^>^^n')v,u' 



K,K 
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where the inertia tensor is defined as 

=^E«A,.<,.'(M«M„').y- (2.13) 



The coupling between the internal and rotational degrees of freedom maintains a complicated structure, i.e. 

B 
2" 

i\/2TB{g{a3d3}'}^ , (2.14) 
where {}^ means the angular momentum coupling to J, M, and the last line corresponds to the A = 3 case. 
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B. Definition of tlie intrinsic system 



We now proceed to the definition of the intrinsic system. We have isolated three rotational variables (the Euler 
angles), but we also have (2A + 1) coordinates a\^i,. It is necessary to reduce their number to (2A — 2). This is 
customary achieved by imposing the condition 



(2.15) 



In fact, the requirement is that the ofF-diagonal elements of J^^^' vanish. These three conditions, which usually are 
non-lineal, in principle reduce by three the number of parameters. The inertia tensor J7„„' can be expressed in terms 
of the coupling to angular momentum two and zero of the collective variables, using the fact that 



J2<lf^l^^'\JM>{MM}i 



(2.16) 



Due to symmetry reasons, the only values of J that survive in the summation related with the inertia tensor are 
two or zero. For A = 3 the explicit expressions for the diagonal components of the inertia tensor are the following: 



Ji=B 



4V7{aa}° + V21 {aa}l - 3Vun{aa}l 



AV7{aa}l 



21 {aa}l + 3Vu^{aa}l 



4\/7 {aa}° - 2V21 {aa}l 



(2.17) 



while the vanishing of the off-diagonal components yield^ 



J7l2 
Jl3 



3SV145R {aa}i 



J23 = 3SV143 {aa}^ = 0. 



(2.18) 



It is convenient to define the auxiliary variables Um — Wo,}m ■ '^^'^ definition of the intrinsic system through the 
cancelation of the off-diagonal components of the inertia tensor is therefore equivalent to the conditions ui — 0, U2 
and u_2 real. It must be noted that similar conditions are held in the quadrupole case. In this case those conditions 
are satisfied if 02,1 = 02,-1 = and a2.2 = 02,-2 and real. In the octupole case, one must satisfy three non-linear 
equations, whose solutions are more complicated. 

The a^iy (in what follows we will concentrate on the A = 3 case, and therefore use a^) must be expressed in terms 
of four independent variables. Calling them X, Y, Z and 7 it is possible to use as a parametrization 

\/3 VS 
03 = [cos7 — sin7]X -I- i[cos7 H — — sin7]y 

02 = — ^ sin7Z 
v2 

oi = -— sin7[X -I- iY] 



ao 



cos "fZ 



(2.19) 



that automatically makes the off-diagonal components of the inertia tensor equal to zero. 

To simplify the notation we define — X,Y, Z and 7k = 7 — f-K for k = 1, 2, 3. In terms of these variables, the 
Um different from zero can be written as 



^ In the quadrupole case J12 is related to 3 {aajj while ^7i3 and are related to {aa}^ . The diagonal components of the 
inertia tensor are related to {aa}[j, {ffla}Q and ^{aa}^- In a similar way in the T = 1 pairing collective description the off 
diagonal components of the inertia tensor can be written as ^12 = — -BSfAA'jj, J7i3 = — -B3?{AA*}^, ^723 = — ^^{AA*}^. 
The diagonal components are related to {AA*}g, 5R{AA*}2 and {AA*}g 
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where V = JiJ2k^I ^osj^. 
If we define 



Uq = T] COS 7 

M2 = ^sin7, (2.20) 



p2 = ^|a,|2 = ^x2[l + 4cos2 7„], (2.21) 



the diagonal components of the inertia tensor can be written as 

J, = b|V_2V2T77C0S7,}. (2.22) 



C. Quantization of the Hamiltonian 

The next problem is that of quantizing the Hamiltonian. There is no unique way to perform this quantification but, 
as in the quadrupole case jl^ , the Pauli prescription can be used. This recipe is designed to give the right answer 
when the generalized variables are transformed to Cartesian coordinates. Given a classical Hamiltonian written in 
terms of variables a^, and their time derivatives, with a kinetic energy given by 

T =-'^G,nn{a)amdn , (2.23) 

mn 

the Pauli prescription replaces the kinetic energy by the operator 

^ = E |Gr^/-|G|^G™"-^ , (2.24) 
2 ^ dam oan 

mn 

where \G\ is the determinant of the inertia matrix Gmn and G™" is its inverse matrix. The volume element to be 
used is given by 

dT=\G\^Wdan. (2.25) 

n 

In our case, the variables are X, Y, Z, 7, (71,(72, 93, and the inertia matrix Gmn is shown in Table ^. 

III. SYMMETRY PROPERTIES 

A. General structure of the wave functions 

The total collective Hamiltonian will contain a potential energy depending on the internal collective variables 
(and its time derivatives) in addition to the kinetic energy terms already discussed. The eigenf unctions of the total 
Hamiltonian can be labeled by the parity, the angular momentum and its projection on the laboratory z-axis (r, / 
and Iz). They can be conveniently expresses as a linear combination 

*J_,JX,y,Z,7,(^,) =^(7^_^(X,r,Z,7)I?,^^_^(c.,), (3.1) 

K 

where the quantum numbers ^, related to the internal variables, remain yet unspecified. This set of eigenfunctions 



constitute a complete set of states and they will be orthogonal using dr specified in Eq. (2.25) as volume element 
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B. Invariants and the potential energy in Bohr's Hamiltonian 



As in the collective descriptions studied before , there are two types of symmetries that arise in the problem 

of collective motion. The first type corresponds to the invariance under rotations. The second type is related to the 
fact that even if the collective variables in the Lab. system are uniquely determined, there are 24 different ways of 
defining a right handed intrinsic system For simplicity, it is commonly assumed that the potential energy does 
not depend on the velocities, and that it should be an analytic function when written in terms of the a^. Naturally, 
it must also be a scalar under rotations, and it is therefore important to search for simple polynomials in terms of 
the that are scalar under rotations. It is well known that in the quadrupole case there are two of the so called 
basic invariants (/?^ and (3'^ cos 87) while in the T = 1 pairing case there are three basic scalars under rotations 
[A^ and A^e^**"^ cos(2r)], but only two if one requires invariance with regard to rotations both in usual space and in 
gauge space [A^ and A*cos^(2r)]. 

In the octupole case it is possible to find in a simple way some of these basic scalars with regard to rotations (we 
use the same notation as Kumar and Baranger The first one is quite trivial (the upper index denote the number 

of bosons needed to construct the scalar) : 



^''^-p' = E"X = -^W}o- (3-2) 



The second and third ones can be constructed in a similar way as in the quadrupole case using the quadrupole 
variables related to the inertia tensor, i.e. 

/(^) ^ {uu}', (3 3) 



J(6) ^ 



|{uu}St|^ = -y^773cos37. (3.4) 



To be sure that one has obtained all the basic scalars it is convenient to count the states with angular momentum 
zero and check that they can be constructed out of the elementary scalars already known. The number of states 
having angular momentum zero for each number of bosons can be constructed in a simple way using the m-scheme. 
Studying the states obtained considering up to 40 bosons it is found that it is necessary to introduce a basic scalar 
formed with ten bosonsQ. Its existence shows how much more complicated is the octupole case as compared with the 
other collective treatments already performed. 

C. Wave function's symmetry properties 

In this subsection we will follow closely Ref. The wave function 'I' must be an analytical function of the 

coordinates a^. The Lab. variables are unambiguously defined while, as Bohr already noted, the intrinsic variables 
am are not, as there are 24 ways of defining a right handed intrinsic system starting from one Lab. system. 

This can be seen more clearly if one defines three basic operators which can be used to transform a given intrinsic 
system in an equivalent one: 

• TZi : rotation through tt around the 1-axis of the intrinsic system. 

TZ2 ■ rotation through ^ around the 3-axis of the intrinsic system. 

TZt, : cyclic permutation of the three intrinsic axes. 



Since TZi =7?,| — TZ'^ — 1, the 24 possible transformations between equivalent intrinsic systems are S{si, S2, S3) = 
7?.^^ 7^.2^72.3^ with < si < 1, < si < 3, and < si < 2. There are then 24 ways of choosing the a„i for a 
given a^, but 4" must be the same for each of the choices, being invariant under the transformation that changes an 
intrinsic system to an equivalent one, i.e. \E' must be invariant under any of the transformations S{si, S2, S3). All 
these transformations can be considered as acting on the coordinates {X, Y, Z, 7, 9, cj), ip). We have therefore that these 
transformations applied to a point 



We would like thank to Prof. J. Blomqvist who called our attention to the existence of this tenth order invariant. 
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. 7^l(X,F,Z,7,^») = iX,Y,Z,j,TT-e,^ + 7r,-ip) 
. 7^2(X, Y, Z, 7, 9, ^) = (X, Y, Z, 7, ^ + f ) 

• Tlj,{X^ Y, Z, 7, 9, (p, ip) — {X, Y, Z, 7, three new Euler angles very complicated) 

These very complicated new Euler angles related to TZ^ correspond to the cyclic permutation of the axis but we 
will not need their explicit form. 

In order to obtain the corresponding transformation Rk of wave functions 5* we must apply the recipe for active 
transformations: 

new wave function at new point = old wave function at old point. 
It is convenient to remember that the application of an active rotation R on a spherical tensor Tim yields 

[Ri9, 0, ij)T]irn = ^™.«(^' (3.5) 

71 

Writing the new wave function at the new point in terms of the old one at the old point one obtains 

R |/, M) = J2 T^kAR^') K) = 5][P^,m(^, 0, n* \I, K) (3.6) 

K K 

a similar relation holds between am and , i.e. am — '/'j ''/')]*Q^^I• 

To obtain the effect of Ri on the 2?^ matrices it is convenient to relate Ri to the operator R'^ that performs the 
same transformation on the Lab. system, i.e. Ri = RR'^R^^, where R corresponds to the rotation {9, </>, ?/>). It is then 
possible to express 

R,Vl,j^i9,(t>,ij) = RR'^R-'v[rj,i9,(t>,^Jj) . (3.7) 



Applying three times Eq. (3.6) on the right hand side one gets 

J2 i?k;^(i?)i?kp(^r')2?p,Q(^"'pLA'(^'^>) = E^k/f(^'<^'^)^kA'(^r'), (3.8) 

KPQ K 

where in the last step the unitarity of the ^ matrices was used. 
Taking into account that the wave functions can be written as 

^}^m{X, Y, Z, 7, 9, 0, V) = E 91k{X, Y, Z, jMi^AO, <P, ^) , (3-9) 

K 

and taking into account that the symmetry of the wave function with regard to Ri , R2 and R3 implies that 

R^^}^MiX, Y, Z, 7, 9, 0, ^) = ^Im{X, Y, Z, 7, 9, V) , (3.10) 
and considering the effect of Ri on the intrinsic and rotational degrees of freedom 

R.Vi,^j,{9A,i') ^T^iLK{&\<l>',^l ^T.^kK'{&,<l>,^^)T^K',K{^^^ (3.11) 

K' 

R^g^^AX, Y, Z, 7) = glAX', Y', Z', 7') , (3.12) 
it is possible to obtain a condition on the part of the wave function that depends on the internal variables: 

gl^,{X,Y,Z,^)^Y.9lK{X',Y',Z',^)V'^,j,{Rr^). (3.13) 



K 



In order to obtain the effect of the symmetry operation on the internal variables, we must first evaluate k{^) 
explicitly. A lengthy yet straightforward calculation yields 
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l^iiM^i) = (-l)'V-if (3.14) 

KLKiR2) - i-ifSM^K (3.15) 

2?W(^3)=2?if,K(f,0,|). (3.16) 

This last matrix is shown in Table |l[ 

From the transformation properties of the coefficients a™, one can deduce the transformation properties for X, Y, Z 
and 7. In Table HI we show all the information related to the transformation properties of the variables, as well as 
the characteristics of the matrices related with I?^ ^, (i?^^^). 

With regard to the spatial inversion V, which is not a symmetry property due to ambiguities in the definition of the 
intrinsic axis, but must be a symmetry of the intrinsic wav e fun ction, one obtains V : {X, Y, Z, 7) —^ (—X, — Y, —Z, 7) 
and j^, = {—1)^Sk.k' that replaces the T)'j^ j^,. in Eq.( 3.13| ). The different transformations will then yield some 
symmetry properties of the wave function: 

i?i : g,V(X,r,Z,7) - {~iyg}^_K{X-Y-Z,j) (3.17) 
i?2 : glAX, Y, Z, 7) = (-*)^".g,V(-F, X, Z, -7) (3.18) 

i?3 : 9J,k{X, Y, 7) - ^ 91k' {X, Y,Z,^+ ^)I?^ ^, (|, |, ^). (3.19) 

K' 

There are some combinations of these operators that yield more useful results, such as 

Rl : g,V(X,y,Z,7) = (-l)^5j^^(-X, -F, Z, 7) (3.20) 
RlR^:glK'{X,Y,Z,-f) = g}^_j,{-X,Y,-Z,j) , (3.21) 

If the wave function is also an eigenfunction of the parity operator with eigenvalue 11, i.e. 

V^}^MiX, Y, Z, 7, 9, V) = Ti^]MX, Y, Z, 7, 9, ^, ^) , (3.22) 

we obtain that 

V : g}AX,Y,Z,^) = \l{-lY g]A-X,-Y,-Z,^) . (3.23) 



IV. THE RIGID ROTOR 



In this section we study the axially symmetric rotor. Ax ial sy mmetry implies that at least two of the diagonal 
components of the inertia tensor must be equal. From Eq. ( 2.22| ) it follows that for 7 = or 7 = -I this condition 



is satisfied. It is also possible to make the three components of the inertia tensor equal by making all the three Xk 
equal (in this particular case 77 = 0). For 7 = one has Ji — J2 while for 7 = f one has J\ = Jz- In both cases the 
remaining component does not vanish. We will study in some detail the 7 = case because it illustrates some of the 
complications and features related to the octupole degrees of freedom. 

We will assume in this section that the 7 degree of freedom is frozen (in particular 7 is equal to zero). It follows also 
from the relation between the collective variables in the intrinsic system (a^) and the intrinsic variables {X^Y^ Z and 
7) given by Eqs. ( ^.19 ) that 7 = implies that the only deformations allowed should have K — Q oy K — ±3, ±6, ... 



It must be noted that when the collective description of the T = 1 pairing excitations was studied in Refs. |gl|,g3[, 
it was found that for studying the rigid rotor it was necessary to retain the gauge angle degree of freedom as the 
rotational Hamiltonian was coupled to it. In this section, even if the rotational Hamiltonian is coupled to the 7 degree 
of freedom we will assume that it is possible to freeze it and we will therefore disregard it completely. In this case the 
two basic quantities and j3 can be written as: 

p2 = 2(^2 + y2) + 5^2^ (4,1) 

/3=iL(-i(X2 + y2) + ^2)^ (4 2) 

and the inertia tensor has a rather simple expression in terms of the intrinsic variables: 
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J3 = 18(^2 + y2) _ (4 3) 

It is convenient to reparametrize the problem in order to introduce p as one of the dynamical variables and use 
"spherical coordinates" by defining two angle variables /i and C 

X = p sin /X cos C/\/2, 
Y — psin/xsinC/-\/2, 

Z = pcos^/V5. (4.4) 
In terms of these variables /? = -^=p'^[—j sin^ + ^ cos^ //]. Now the Gmn inertia matrix will be a six by six matrix 



as shown in Table IV 



We will assume the adiabatic hypothesis for p and p. But it must be noted that, as C is a cyclic coordinate, we 
cannot use for it the adiabatic hypothesis and we must, therefore, consider the coupling of C with qi, 52 and 53 in an 
exact way. For this purpose we use the conjugate momentum method, i.e., we introduce 

dT 3 

Qi = — =6p2[l- sinVlgi (4.5) 
oqi 4 

dT S 

Q2 = ^=6p2[l--sinV]'Z2 (4.6) 

dT 1 ■ 

Q3 = ^ = sin^ ^[^3 - 3C] = ^ (4.7) 

dT 1 ■ d 

PC = -^^ ~3p2 sin2 p[q, - -C] = -I— = . (4.8) 



It must b e no ted that Qi satisfy the commutation relations given by Eqs.( |2.8| ). Besides, the relation between Eq.(4.7) 
and Eq. ( |4.8| ) provides the auxiliary condition: 



iv: = -3to(^ = 0,±3,±6, ... (4.9) 
The kinetic energy can be written in terms of these new variables (3„ as 

^-*=E^Q'' (4-10) 
where the components of the moment of inertia are given by 

Ji = :72 = ^P^[l + 3cosV] 

J-i = ^p^ sin^ p . (4.11) 
The total rotational energy can be written as 

The structure of the full "rotational band" is now more complicated than for the quadrupole case. The part of the 
band that has K = has a structure similar to that of the usual rigid quadrupole band, but the part of the band 
related to the states with K ^ has energies determined also by p and p that are related with the energy of the 
first member of the ground state rotational band. Note that p represents deviations from perfectly axisymmetric 
deformations, i.e. ioi p — the deformations have jTs = and therefore K = 0, only for p ^ can nonzero K be 
found. For p = the rotational band has the usual structure (all states have K = 0, as states with K have an 
infinite energy) . When the states with angular momentum equal or larger than 3 are considered it is found that there 
are more than one state belonging to the ground state band with this angular momentum: one with K ^ and parity 



(— 1)^, and two with K = 3 and energies given by Eq. ( 4.12 ). Each time that / = \K\ = 3n (with n integer) two new 



states appear (i.e. we will have five states with 8 > / > 6, seven states with 11 > / > 9, etc.). Figure ^ shows the 
lowest energy states for if 0, 3, 6 as a function of p. 
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FIG. 1. Spectrum of a symmetric (7 = 0) rigid rotor for /sT = 0, 3, 6 as a function of the deformation parameter /u. 




V. SUMMARY 



In the present paper we have obtained a parametrization of the octupole collective variables that guarantees that 
the inertia tensor is diagonal. The relation between the components of the octupole intrinsic variables and those 
defined parameters is given by the equations 

as = [cos 7 — ^ sin7]X + i[cos7 + ^ sin 7]!^ 

02 = -j^sin'yZ 

ai = ^ sin^lX + iY] 
ao = ^cosjZ 

In addition, we have identified the way in which the intrinsic variables (xfc and 7) transform with regard to the 
symmetry operations Ri. These can be summarized as follows: 

a) The variables Xk transform in a similar way as the coordinates of a vector related to the intrinsic axis, 

b) 7 transforms in the same way as in the quadrupole case. This is not at all strange due to the particular relation 
of 7 to the inertia tensor. 

The first condition means that the parameters X, Y and Z can be considered as positive, as all the other signs can 
be obtained by relabeling the axis. The second statement, as in the quadrupole case, tell us that all possible shapes 
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are contained in the range < 7 < ^. 

A general octupole field can be written in an alternative way ||l6| as V3 = X^m^m^Sm = ^0^2 + J2'i=i + 
X]i=i ^2ii)F2{i)i where Ai, Fi{k) and F2{k) are related to the irreducible representations of the octahedron group. In 
Ref. ||l^ it was proposed that an appropriate parametrization can be obtained imposing the cancelation of e2(*)Vi. 
Using the parametrization that we are proposing, these parameters have now a rather simple structure: 

£0 = 

ei(l) = -sin(7)Z 
6i(2) = -sin(7-f )X 
ei(3) = sin {l~^)Y 
£2(1) = \/5cos(7)Z 
62(2) = x/5cos (7- 
£2(2) = -V5cos {l^)Y . 

It must be noted that, for the octupole degrees of freedom, the existence of axial symmetry does not necessarily 
imply that the rotational band must have K = 0. The study of the rigid rotor has illustrated the richness of the 
octupole rotational bands. 

Last but not least, it is worthwhile noticing that if one takes into account simultaneously the quadrupole and 
octupole degrees of freedom, and for each one of them the condition that the inertia tensor has to be diagonal is 
imposed, the intrinsic system for both degrees of freedom will not be the same: each degree of freedom will have its 
own "intrinsic system". It must be remembered that the Euler angles associated to each intrinsic system are just a 
way of labeling three of the five (in the quadrupole case) or seven (in the octupole one) collective dynamical variables. 

Discussions with R. P. J. Perazzo and J. Fernandez Niello are gratefully acknowledged. This work has been supported 
in part by the Carrera del Investigador Cientffico y Tecnico, by PID iV°4547/96 of the CONICET, Argentina, PMT- 
PICT1855 of ANPCYT and Grant Ex-055 from the University of Buenos Aires. CW was supported in part by the 
University of Buenos Aires. 
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X 


Y 


Z 


7 


qs 


q2 


qi 


X 


l+4cos^7i 








-2Xsin(27i) 


4Y(l+2cos7i cos 72) 


-4Z(l+2cos7i cos 73) 





Y 





l+4cos^72 





-2Ysin(272) 


-4Z(l+2cos7i cos 72) 





4Z(H-2cos72 cos 73) 


Z 








l+4cos^73 


-2Zsin(273) 





+4X(l+2cos7i cos 73) 


-4Y(l+2cos7i cos 73) 


7 


-2Xsin(27i) 


-2Ysin(272) 


-2Zsin(273) 




-6V3xy 




-6V3yz 


qs 


4Y(l+2cos7i cos 72) 


-4Z(l+2cos7i cos 72) 





-&^ZXY 










q2 


-4Z(l+2cos7i cos 73) 





4X(l+2cos7i cos 73) 







J2 





qi 





4Z(l+2cos72 cos 73) 


-4Y(l+2cos72 cos 73) 


-6V3yz 








Ji 



TABLE II. Matrix elements of 2?M,K(f , 0, f ) 



M\K 


3 


2 


1 





-1 


-2 


-3 


3 


i 
8 


8 




"^"^^^ 


,• Vl5 

8 


8 


i 
8 


2 


v& 

8 


4 
8 


Vio 

8 





Vio 

8 


4 


V6 


1 


»• VIS 
?- 


VIO 
* 8 


^ 


• \/l2 
' 8 


i 


vio 

* 8 


»• V15 





8 





^12 
8 





V^2 

8 





V20 

8 


-1 


^ 


»• Vio 

' 8 


i 


,•^12 
* 8 


i 
8 


»• Vio 

' 8 


Vio 

gy- 


-2 


8 


4 







\/l(J 




4 


V6 
8 


-3 


i 
8 


^ 

' 8 


'~8~ 




'~8~ 





i 
8 



TABLE III. Action of the different transformation operators on the corresponding intrinsic variables. Also shown are the 
explicit representation of 'Dk^k'I^T^) 



Variable\Trans formation 


Ri 


ii2 




X' 


X 


-Y 


Z 


Y' 


-Y 


X 


X 


Z' 


-z 


Z 


Y 


7 


7 


-7 


7+^ 


9 


TT-e 


6» 


Complicated 




TT + 


</> 


Complicated 








Complicated 








^K,K' V 2 ' 2 ' ^l' 



TABLE IV. Inertia matrix Gmn for the 7 = case 





p 




c 


91 


92 


93 


p 


1 

















A 



















C 








siv? fi 


o 2 • 2 

—Sp sin p 








91 








—Sp sm fx 


9p^ sin'' /Lt 








92 














6p"[l-f sinV] 





93 

















6p^[l-|sin>] 
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